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We study a Hartree ensemble approximation for real-time dynamics in the toy model of 1+1 dimensional scalar 
field theory. Damping behavior seen in numerical simulations is compared with analytical predictions based on 
perturbation theory in the original (non-Hartree-approximated) model. 



1. Introduction 

Dynamics of quantum fields in real time is a lot 
more complicated than statics in imaginary time 
and one has to make approximations before giving 
the problem to the computer Q . As an improve- 
ment on classical dynamics we are studying the 
Hartree approximation. The Hartree equations of 
motion for the ip 4 model in 1+1 dimensions are 
given by 

Cp x = Aip x - (n 2 + Xipl + 3AC X ) tp x , 
fS = A/«-( M 2 + 3A^ + 3AC a )/«, 
C x = £[(l + 2n ")|/«| 2 ], 

a 

where tp x = {<p x ) is the mean held and the f x 
are a complete set of mode functions used for the 
parametrisation of the operator field (p x in the 
Hartree approximation, 

a. 

The initial conditions are specified by a suitable 
choice for ip x , ip x and the mode functions, as well 
as the initial occupation numbers rig = {a) a a a ). 

Hartree-like approximations have been widely 
used, but when the mean held is homogeneous 
they do not lead to proper thermalization. This 
may be ascribed to lack of sufficient scattering in 
the interaction term involving a spatially constant 
C x . We try to improve on this Q by writing a ho- 
mogeneous initial density matrix as a superposi- 
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Figure 1 . Numerically computed auto-correlation 
functions ln|.F m f(i)| versus time mt, with m the 
temperature dependent mass. The coupling is 
weak, A/to 2 = 0.11 and the temperature T/m w 
1.4 for the smaller volume (with significant devi- 
ations from the Bose-Einstein distribution) and 
Ri 1.6 for the larger volume (reasonably BE). 

tion of gaussian pure states. The Hartree approx- 
imation is then applied to realizations of this en- 
semble, which are typically inhomogeneous. This 
allows for scattering of the modes ('particles') via 
the fluctuating mean field and redistribution over 
various momenta. To test for thermalization we 
start out of equilibrium with 'the modes in their 
vacuum' (tiq = 0) and compute equal-time- and 
auto-correlation functions. For the resulting par- 
ticle distributions, see Here we concentrate 
on damping phenomena. 
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We found funny modulations in auto- 
correlation functions. Fig. |l| shows such mod- 
ulations on top of a roughly exponential decay. 
The correlation function corresponds to the zero 
momentum mode of the mean field, 



F iai (t) = j dx <p(x,t)tp(0,0) conn , 

where the over-bar denotes a time average, taken 
after waiting a long time for the system to be in 
approximate equilibrium. As discussed in [^| , this 
equilibrium is approximately thermal. The func- 
tion F m f(t) is analogous to the symmetric corre- 
lation function of the quantum field theory: 

F(t) = J dx (\{$(x,t), 0(0,0)}) conn- (1) 

Does F(t) also have the modulations? 

2. Calculation of F(t) 

The function F(t) can be expressed in terms of 
the (zero momentum) spectral function pip ), 



F(t) = f 

J —i 



dp 
2tt 



-i P u t 



1 1 



e P°/T_i 2 



P(P°)(2) 



which is determined by the retarded selfenergy 

£(p°), 



p(p°) = 



-2ImE(p°) 



,2 _ „2 



p2 + ReI](pO)]2 + [l m S(pO)] : 



(3) 



The selfenergy can be calculated in perturba- 
tion theory. The relevant diagrams are shown in 
Fig. H. The one-loop diagram is present only in 
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Figure 2. Diagrams leading to thermal damping. 



the 'phase of broken symmetry'. It leads to damp- 
ing for frequencies p\ > 4m 2 , which are irrelevant 
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Figure 3. Real (left) and imaginary (right) part of 
— i In (i/z + \J\ — l/z 2 ^j = arcsin(l/z), for z = 
(p° +ie- m)/ (0.13m) and e/m = 1(T 4 . 



for the quasiparticle damping at y\ 



(There 



is really only a symmetric phase in 1+1 dimen- 
sions, but this is due to symmetry restoration 
by nonperturbative effects which will not oblit- 
erate the damping mechanism described by the 
diagrams.) So we concentrate on the two-loop di- 
agram. This is given by the sum of two terms, 
Ei + E2. The first has an imaginary part corre- 
sponding to 1 <-* 3 processes, requiring p 2 , > 9m 2 , 
so it does not contribute to plasmon damping. 
The second is given by 



9A 2 



dp 2 dp 3 



16tt 2 J E1E2E3 
(1 + ni)n 2 n 3 - m(l + n 2 )(l + ^3) 



p u + 1 
[(P° + ie) 



E\ — E% - 
-(p°+ie) 



E% 



(4) 



where A is the coupling constant (introduced as 
£1 = -A^ 4 /4 ), and 



Ei = \Jm? + (p 2 +pz) 2 , E 2 ,z = \Jm 2 + p\ z , 
n, = [exp(E t /T)-l}-\ 1 = 1,2,3. 

Its imaginary part corresponds to 2 <-> 2 pro- 
cesses, which contribute to plasmon damping (the 
regions near po = ±m ) . 

Now the formula for the thermal plasmon 
damping rate (at zero momentum) in terms of 
the retarded selfenergy, 



7 = — ImE(m)/2m, 



(5) 



leads to a logarithmically divergent answer. This 
is a collinear divergence which is absent in more 



Figure 4. Plot of Re £(p )/(9A 2 /16tt 2 ) obtained 
by linear extrapolation e = 0.02, 0.01 to zero, 
together with a matching to the logarithmic sin- 
gularity [T = m = 1). 



than one space dimension. Inspection shows that 
the singular part of £2 is given by the non- 
relativistic region of the integral in (^). Using 
polar coordinates P2 = pcoscf), p% — psin</> this 
non-relativistic (jp < n <C to) contribution is pro- 
portional to 

P d P / — 

Jo 




le - 



k 2 /2m 



where p° w to. This function is plotted in Fig. [| 

A natural way out of the difficulty posed by the 
divergent £(to) may be to continue the selfenergy 
analytically into the lower half of its second Rie- 
mann sheet, p° — > to — ij, and replace (^|) by the 
improved definition 

to 2 — (m — i"f) 2 + £(m — £7) = 0. (7) 

For weak coupling A/m 2 C 1 we then get the 
equation 

7 9A 2 e m ' T 



16nm 4 



( e m/T _ X y 



In- 
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c(T) 



(8) 



Figure 5. As in Fig. 4 for the imaginary part. 



where the constant c has to be determined by 
matching a numerical evaluation of E to the form 
(||) for p° w rn. 

We evaluated £2 in (§) for T = to by numerical 
integration with e/m = 0.02, 0.01 and linear ex- 
trapolation e — > 0. The result is shown in Figs. 4 
and 5, together with a matching to the logarith- 
mic singularity, giving c ss —0.51. For example, 
Eq. (||) now gives 7/771 = 0.061, for A/to 2 = 0.4. 

To see how well this 7 describes the decay of the 
correlator F(t) we evaluated this function directly 
from (|^) and (|^). The divergence in Im£(p°) at 
p° = to leads to a zero in the spectral function 
p(p°). So is there a peak at all in p(p°)? Fig. 6 
shows what happens: the 'usual' peak has sepa- 
rated into two twins! 

Fig. 7 shows the resulting F(t). The effect of 
the double peak is indeed an oscillating modula- 
tion on top of the roughly exponential decay. The 
decay corresponding to exp(— 7*), with 7 given by 
(||), is also indicated in the plot: it does not do a 
good job in describing the average decay beyond 
the first interference minimum. The 'Twin Peaks' 
phenomenon implies that the usual definition of 
damping rate (0) is unreliable in 1+1 dimensions. 

Fig. 8 shows the result of a calculation of F(t) 
with parameters taken from the numerical simu- 
lation in Fig. 1 with the larger volume. In this 
case e was kept finite, e/m — 0.005, which may 
be more realistic since one expects anyway the 
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Figure 6. The spectral function p(jP) near p° = 
777. = 1 corresponding to the selfenergy shown in 
Figs. 4, 5 (T = m, A = 0.4m 2 ). 



infinity in Im £ to be smeared out by damping 
effects in the propagators in Fig. g. Figs. 1 and 
8 are reasonably similar, but how to compare the 
the average slopes is somewhat ambiguous. 

3. Outlook 

Summarizing, we are encouraged by the simi- 
larities in the qualitative features of the numerical 
and analytical auto-correlation functions. Quan- 
titatively, the damping times are also of the same 
order of magnitude (105 vs 67 m _1 ), but we did 
not really accurately compute the relevant auto- 
correlation function yet in our simulations, which 
will require a lot more numerical effort. 
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Figure 7. Plot of In \F(t)\ versus rat for T = m, 
A = 0.4m 2 . 




Figure 8. ln|F(i)| versus mt with A/777 2 = 0.11, 
T/m = 1.63, corresponding to Fig. 1. The line is 
given by exp(0.55 — t/67). 



